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i-Q 1. Introduction 

a The structure for spatial permutations consists of a large box A C M'', a large 

, , number N of points in A, and permutations of those points such that all permutation 

jumps remain small. The relevant parameter is the density p = N/\A\. In many 
CO models there is a critical density Pc that corresponds to a transition from a phase 

^ with only finite cycles (when p ^ p,,) to a phase where a nonzero fraction of points 

^I belong to infinite cycles (when p > pc). The goal of the present article is twofold. 

*vj First, we prove that such a transition occurs in a class of models of spatial random 

psj permutations with cycle weights. Second, we show that the cycle structure of 

I I infinite cycles satisfies a Poisson-Dirichlet law. 

f^ The main motivation for our models comes from the interacting Bose gas of 

f^ quantum statistical mechanics. The possible relevance of long permutation cycles to 

'"^ Bose- Einstein condensation was pointed out by Matsubara [T7] and Feynman [TU] . 

^ Siito made important clarifications for the ideal Bose gas, showing in particular 

that long cycles are macroscopic [HI [20]. It is a notoriously difficult problem 

to prove Bose-Einstein condensation. Another problem, that is related but not 
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^ subordinated, is to understand how the critical temperature is modified by particle 

interactions. In the recent article [5], we derived (non- rigorously) a model of spatial 
permutations where the original interactions between quantum particles have been 
replaced by cycle weights. The simplified model retains some features of the original 
model, as they have the same free energy to lowest order in the scattering length 



of the interaction potential. We then used the formula (2.9) below for the critical 
density. The validity of this formula for the model with cycle weights is proved in 
the present article. 
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2 VOLKER BETZ AND DANIEL UELTSCHI 

Models of spatial permutations are also attractive per se. They have both specific 
and general features. One general feature that is especially striking is the Poisson- 
Dirichlet law for the distribution of cycle lengths. The literature on the subject is 
huge, see e.g. [TJ[T31[Ilj for a sample. The Poisson-Dirichlet distribution is expected 
to make an appearance in other models with spatial structure and permutations 
such as the random stirring model jl3|, I21j . This was proved recently by Schramm 
on the complete graph |18j; see also Berestycki \2, for several useful observations 
and clarifications. 

The models considered here are "annealed" in the sense that spatial positions 
vary and they are integrated upon. Annealed models are both simpler and more 
relevant for the Bosc gas. But the "quenched" models, where the positions are fixed 
and chosen according to a suitable point process, look very interesting in probability 
theory. One conjectures that long cycles satisfy the same Poisson-Dirichlet law 
as in the annealed model — the only difference being the critical density. This 
is supported by numerical evidence [111 I15| . An unrelated but very interesting 
problem is the complete description of Gibbs states, involving crossing fluxes that 
depend on the boundary conditions. Such a description has been recently achieved 
by Biskup and Richthammer in the one-dimensional model [7]. 

2. Setting & results 

The state space of the (annealed) model of spatial permutations with cycle 
weights is ^a,n — A^ x Sn, where A C K'' is a cubic box of size L, and Sn 
is the symmetric group of permutations of N elements. We denote by | A| = L'' the 
volume of A. We equip 51a, w with the product of the Borel cr-algebra on A^ and 
the discrete cr-algebra on Sn- We introduce a "Hamiltonian" and its corresponding 
Gibbs state. Namely, the Hamiltonian is a function H : ^Ia^n — ^ K U {oo} that we 
suppose of the form 

N 

H{x, tt) = ^ £_{xi - x^(i-j) + ^ airi{n). (2.1) 

i=i e^ 1 

Here, x = {xi, . . . , xn) G A^ and tt G Sn- We always suppose that e~^ is contin- 
uous with positive Fourier transform, and that it is normalized, /g^ e~^'^-* dx — 1. 
Notice that ^ is allowed to take the value -|-oo, and that positivity of the Fourier 
transform implies that ^{x) = ^{—x). The cycle weights ai,Q;2, . . . are fixed pa- 
rameters. Finally, reij:) denotes the number of ^-cycles in the permutation n. 

Boundary conditions are not expected to play a prominent role here, and we 
therefore choose those that make proofs simpler. These are the "periodized" bound- 
ary conditions, where we replace ^ by the function ^a, defined by 



zeV 



4^1 imnlips that p ?a(2:) 



IS 



The normalization assumption /» e ^'^ = L^ e ^ = 1 implies that e 

finite for at least almost every x- When e~'> has bounded support with diameter 

smaller than L we recover the usual periodic boundary conditions. We let Ha be 



as H in (2.1 1, but with ^a instead of ^. The Gibbs state is given by the probability 
measure 

Prob(da;, tt) = -^ e^^'^^^^'^) dx (2.3) 
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on flA,N, where dec is the Lebesgue measure on A^ and F is a suitable normaUza- 
tion, namely 

r = ^ y / e-"'^^''''''^ dx. (2.4) 

In typical realizations of the system, points are spread all over the space because 
of the Lebesgue measure that prevents accumulations. The lengths of permuta- 
tion jumps \\xi — a;7r(i)|| stay bounded uniformly in A because of the jump weights 
Q-iixi-x^^i-,) ^ fpj^g lengths of permutation cycles depend on the density of the sys- 
tem. For small density, points are far apart and jumps are unlikely, which typically 
results in small cycles. But as the density increases, points have more and more 
possibilities to hop, and a phase transition takes place where "infinite" cycles ap- 
pear. The cycle weights modify the critical density and also the distribution of 
cycle lengths, see below. The model is illustrated in Fig. [TJ 




Figure 1. A typical realization of a spatial permutation. As 
|A|,A'^ -> cxD, the jumps remain finite but the cycle lengths may 
diverge. 

This model arises naturally from the Feynman-Kac representation of the dilute 
Bosegas. The jump function is then ^(x) = ^|la;|p (plus a normalization constant), 
with /3 the inverse temperature of the system. Notice that if the original quantum 
system has periodic boundary conditions, we get the periodized Gaussian function. 
Cycle weights were introduced in [3] as a crude way to account for the particle 
interactions. But the calculations of [S] suggest that the cycle weights can be 
chosen so that the model describes the Bose gas exactly in the dilute regime. We 
do not write here the precise formula for the weights, but we observe that they 
satisfy the asymptotic aj = —a{l — 0{j^^'^)), so that aj converges as j — >■ cx) fast 
enough for our purpose. 

We are solely interested in properties of permutations and we introduce random 
variables that are functions on Sn rather than flA,N- Let i'^^^n) , fS'^'> (tt) , . . . denote 
the cycle lengths in non-increasing order, repeated with multiplicities. We will prove 
that, above the critical density, the cycle lengths scale like N and they converge 
in distribution to Poisson-Dirichlet. The latter is conveniently defined using the 
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Grijfiths-Engen-McCloskey distribution GEM(0), which is the distribution for 

(Xi, (l-Xi)X2,(l-Xi)(l-X2)X3, ...), 

where Xi, X2, ■ ■ ■ are i.i.d. beta random variables with parameter (1,^); that is, 
Piob{Xi > s) — {1 — s)^ for ^ s ^ 1. The Poisson-Dirichlet distribution PD(0) is 
the law obtained by rearranging those numbers in non-increasing order. See [U [14] 
for more information and background. In the sequel, we say that a sequence of 
random variables Yn , !« , ■ • • converges in distribution to Poisson-Dirichlet as 
n — )■ 00 if, for each fixed k, the joint distribution oiYn , . . . , Yn converges weakly 
to the joint distribution of the first k random variables in Poisson-Dirichlet. This 
is denoted 

(rW,y,p),...)^PD(0). (2.5) 

As already mentioned, we make the important assumption that the jump func- 
tion has nonncgative Fourier transform. This allows to define the "dispersion rela- 
tion" e{k), k E M'*, by the equation 

e-"^''') = / g-27rife.x ^-^iix) j^^ ^2.6) 

Notice that e{k) is real, e(0) = 0, and e{k) > for all fc 7^ 0, and limj|;;j|^oo £{k) = 00 
(by Riemann-Lebesgue) . In order to avoid pathological cases we assume that e"'^*^'^-' 
is uniformly continuous on R''. We also suppose that e{k) ^ a||A:|j'' for small k, for 
some a > and rj < d. It is easy to see that £{k) is always greater than a||fc|p for 
small k, so the latter assumption always holds in dimensions d > 2. Among possible 
jump functions other than Gaussians, let us mention e"'^'^-' = const (|a;|-|-l)~''' with 
l<7<2ind=l, for which 77 = 7 — 1. As for the cycle weights, we consider three 
cases: 

(i) limj_j.oo Oj — a with a > 0, and ^ \aj ~ a\ < cx). 
(ii) limj^oo ctj = a with a ^ 0, and ^ ^\aj — q| < 00. 
(iii) aj — jlogj with 7 > 0. 

We now introduce the fraction v of points in infinite cycles. It is obvious that 
finite systems can only host finite cycles, so the definition of i' must involve the 
thermodynamic limit. Given a finite number K, let vx denote the fraction of points 
in cycles of length larger than K. Precisely, 

i^K= limini e(^^ y ^(*'V (2.7) 

Here and in the sequel, the limit |A|, A^ — ^ cx) means that both go to infinity while 
keeping the density p = A^/|A| fixed. This is the standard thermodynamic limit of 
statistical mechanics. We then define 

J/ ^ lim i^K- (2.8) 

This limit exists since (i'k) is decreasing and bounded. Let D^ denote the lim sup 



of (2.7). We expect that Pk = vk but wc do not prove it. On the other hand, we 
will prove in Section [5] that vk also converges to v as K -^ 00. 
Next we introduce the critical density by 

e-^'^^'^) dk. (2.9) 



j> 1 
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It follows from our assumptions that the critical density is finite. Indeed, the 
numbers e~"j are bounded, so Pc is bounded by the integral of a geometric series, 
/ ce(fc^)-i ' which is finite. 

We propose now two theorems that confirm that pc is indeed the critical den- 



sity of the model, at least in several interesting situations. The formula (2.91 is 
presumably valid beyond the cases treated in this article, but the precise extent 
of its validity is not clear. The first theorem states that macroscopic cycles occur 
precisely above the critical density, and that they obey the Poisson-Dirichlet law. 

Theorem 2.1. Assume that aj -^ a as described above. Then 

(a) the fraction of points in infinite cycles is given by 

V — max! 0, 1 ^ j ; 

(b) when v > Q, i.e. when p > p^, the cycle structure converges in distribution 
to Poisson-Dirichlet: As |A|, A^ -^> cx) we have 






PD(e""). 



Such a law was already observed in absence of spatial structure, and when the 
cycle weights are constant. This case is known as the Ewens distribution, see e.g. 
[21 [121 H] ■ Results about weights that are asymptotically Ewens can be found in 
[161 [^. Spatial permutations with small cycle weights, i.e. when the limit is a = 0, 
were studied in [4 . 

The second theorem concerns cycle weights that diverge logarithmically — it is 
somehow the limit a — >■ oo of Theorem [2.1| Cycle weights have a striking effect as 
a single giant cycle occurs above the critical density! This is in accordance with a 
similar observation for non-spatial permutations [Sj. 

Theorem 2.2. Assume that aj = j\ogj with 7 > 0. Then 

(a) the fraction of points in infinite cycles is given by 

V = max U, 1 

^ P 

(b) when v > Q, i.e. when p > p^, there is a single giant cycle that contains 
almost all points in infinite cycles: As |A|, A^ — > 00 we have 

vN 

The rest of this article is devoted to the proof of the results above. We reformu- 
late the problem in the Fourier space in Section Isl following Siito [20] . The model 
involves a measure on occupation numbers of Fourier modes, and of random per- 
mutations of those numbers. In Section [4] we obtain information about occupation 
numbers using techniques of Buffet and Pule [5] , and using certain estimates of our 
recent joint work with Velenik [6]. Random permutations within each mode involve 
the cycle weights and are thus similar to those studied in [6 . Combining all those 
results allows us to prove Theorems |2.1| and |2.2| in Section [Sj 
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3. Random permutations and Fourier modes 

The goal of this section is to introduce an alternative model of random permuta- 
tions that involves Fourier modes, and that has the same marginal distribution on 
cycle lengths. Let A* ~ jTJ^ be the space dual to A in the sense of Fourier theory. 

3.1. The marginal distribution of cycle lengths. Recall that the cycle struc- 
ture of a permutation tt e Sj^ is the sequence of cycle lengths £ — {i^^\£'^^\ . . . , i?(™)), 
with ^(*^ ^ £(«+!) and l^""^^ ^ 1; the number of cycles m depends on tt, 1 ^ m ^ A^. 
Those numbers form a partition of {1, . . . , A^}. Another way to write £ is to in- 
troduce the occupation numbers r = (ri, . . . jT^v), where rj = #{i : £'*•' = j}. We 
always have 

m N 

J2i(^)=Y.jr^^N. (3.1) 

One should not confuse the occupation numbers r with the occupation numbers 
n = {iT-k) to be introduced later; they are not related in any direct way. 



Proposition 3.1. The marginal of the probability measure (2.3) on occupation 
numbers is 



N 

1 XT- 1 / e""^ 



P-Mr-)4n;:^(V^ 



e-js(fe) 



-1. J. J* . I \ a 

j = l ^' -^ fcsA* 



with Y the normalization of ( |2.4[ ). 

Proof. The marginal probability on permutations is 



1 r « (3-2) 



N\Y ja« 

We observe that integrals factorize according to permutation cycles. The contribu- 
tion of a cycle of length j is (with yj+i = j/i) 

e-"^ f e-^-i«*(«'-2''+i)dyi...dyj= e-"^|A| ^(e-«)*''(Lz). (3.3) 

To see the equality in (3.3 1, we start with the right hand side. Using the defini- 
tion of the convolution, writing |A| = J.dyi, and shifting all the variables in the 
convolution integrals by j/i gives 



= y" / dj/i / dy2 ■ • • dyj e 



-i(Lz-y2+vi) g-{(i;2-J/3) Q~iiVj~yi) 



= y" / dyi y f dy2.-.d%e-«(^i-^^+-^(^i-^^» e~«(2^^"^^+^(^^-^^)) 
zisz" 22,...,ZjeZ'' 

(3.4) 

The last equality is obtained by decomposing the domain of integration M'^ into 
cubes A + Lz with z G Z'^ and then changing variables in the integrals so that all 
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the boxes become centered at 0. We now change to summation index: Zj = Zj, and 
Zi — Zi — Zij^i for i < j. It is easy to see that this is indeed a bijection on {Z'^y . 
Summing over Zi instead of Zi now gives the left hand side of (3.3 1. 

The Fourier transform of (e~'^ )*■' is e"-''^'^'^^ . The Poisson summation formula 
states that 

E /(^^) = ^ E /(fc)' (3-5) 

where / is the Fourier transform of /, whose precise definition can be found in Eq. 
([2^1). We then get 



1 ^ 



feSA* 



je{k) 



Tj (tt) 



(3.6) 



All permutations of a given cycle structure have the same probability, and there 
are 

(3.7) 



n,/ 






elements in the cycle structure defined by r. We get the claim by multiplying the 
above probability by this number. D 

3.2. Decomposition of permutations according to Fourier modes. We de- 
note by Ti = (rifc) the occupation numbers indexed by /c G A*, and by AfA,N the set 
of occupation numbers such that X^fepA* '^fe ~ ^- Next, we introduce permutations 
that are also indexed by Fourier modes, tt = (tt^.). Let A4a,n be the set of pairs 
(n, 7r) where n G N'a.n and tt = (tt^) with tt/j e 5„^, for each fc G A*. We introduce 
a probability measure on Ma.n'- 



Prob(n) = ^ n ^"'^'^"' ^"- 



feeA* 



with 



^" = ;;i E 



-^r. 



i''j('^) 



(3.8) 



(3.9) 



7re5„ 



and /lo = 1. We will check later that the normalization Y is the same as given in 
(2.4). Then we introduce the probability of a pair (n, tt) by 



Prob(n,7r) = ^ [] - 

y J- J- rii. 



Y -- nfe! 

keA' 



^{k)nk-J2j J laj'-jC'Tfc) 



(3.10) 



Notice that (3.81 is the marginal of (3.101 with respect to tt. The conditional 
probability Prob(7r|n), where TTk G 5„^ for all fc, is given by 



Prob(7r|n) = ]J 



keA' 



nk^-Kk 



-.-Ej > iQjrj(7rfc) 



(3.11) 



That is, given n, each Hk is independent and distributed as nonspatial random 
permutations with cycle weights (see Eq. (5.1 ) below). Given tt, let Tj — ^j, rj{'Kk). 



Proposition 3.2. The marginal of the probability measure (3.10) with respect to r 
is identical to the marginal of the probability measure (2.3). 
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Proof. We check that the marginal of (3.10) gives the formula of Proposition 3.1 
For this, let r be a collection of occupation numbers, and write (rj^) : r for the set 
of all integers rjk {j > 1, /c G A*) such that J^k ^jfc ~ ^j ^'^^ ^^^ .?• Then, 



P..„bw4i: i: n(ii 



-ir{k)nk-J2i ajrj{TTk) 



{rjk)-r (n.,7r): fceA* 
rj{-!Tk)=rjk 






(3.12) 



{r.ik)-rkeA 



We have summed over jTk that are compatible with r^k, using the formula (3.7 1 for 
the number of elements. The bracket above factorizes according to j. Using 



n 

/cGA* 



/ e~"3 Y' 



(3.13) 



and writing, for fixed j ^ 1, (rjk) '■ tj for the set of integers Vjk with ^j, rjk 
we get 



Prob(r) = T7 n 



o> 1 



J 



" E n ; 



^-je(fc)rjfc 



(rjfc):r3 fceA* 



''ifc! 



For each fixed j, the multinomial theorem gives 

Q-je(k)rjk I 



E n 



r-jfc 



(E 

fceA' 



-j£(fe) 



(3.14) 



(3.15) 



Then Prob(r) is indeed given by the formula of Proposition 3.1 This also proves 
that Y is the correct normalization that makes (3.8 1 and (3.10) probability mea- 
sures. D 

4. Properties of occupation numbers 

We study in this section the probability measure of occupation numbers of 
Fourier modes, Prob(n), that is defined in (3.8). We show that the typical n 
has the following properties: 



"0 _ 

N ~ 



max(0,l-£f); 



• ~k X]o<llfcll<i5 ""fc ^^ small when 5 is small. 

• For all 5 > 0, ;^ X]||A;|| > 5 nklni,>M is small when M is large. 



The behavior of the normalizations /i„ defined in (3.9) play an important role. 

;ants 

(4.1) 



We assume that ft,„ grows or vanishes at most polynomially, i.e., there are constants 
C and K such that for all n ^ 1, 



We also need that certain ratios of /i„ be bounded. Precisely, for s > 1, let 



C{s) 



sup — 

771, n ^ 1 '^n 
n/s<m<sn 



(4.2) 



We assume that C{s) is finite for any s. 

Those properties have been verified in [B] when aj — ?> a and aj = jlogj. Indeed, 
one finds /i„ ^ n^'', with r = e^" — 1 in the first case and r = — 1 — 7 in the second 
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case. The results of the present article actually extend to other cycle weights, as 
long as Eqs (4.1 1 and (4.2) hold true. 

The radius of convergence of the generating function of ft,„ is equal to 1. We 
have the following identity for all 7 > 0: 

J2 e-^" K = exp ^ i e-^^-"^ . (4.3) 

n ^ j" ^ 1 

See [H [TB] for the proof. 

4.1. Macroscopic occupation of the zero mode. We use a strategy that is 
inspired by Buffet and Pule in their study of the ideal Bose gas [8]. It consists 
in looking at the Laplace transform of the distribution of ^. Let Y{N) be the 
normalization of Eq. (2.4 1. We now put the explicit dependence on N because it 
is going to vary. Notice that Y{N) also depends on A, but the domain is fixed 
throughout. 
We have 

n.o=3 

with 

^W= E n e-^e^^-^/^n.. (4.5) 

nSA/A.jv feSA* 
no=0 

Notice the relation 

N 

Y{N)^Y.^,Y{N~]). (4.6) 

3=0 
We are often going to interchange infinite sums and products. Let A^ be the set 
of finite sequences of integers. Notice that M is countably infinite. The following 
lemma is easy to prove and sufficient for our purpose. 

Lemma 4.1. Let a{k, n) be a nonnegative function such that a{k, 0) = 1 for all k. 
Then 

E n "(^'"fe) = n ( E °''^^''' 

(It is possible that both sides are infinite.) 

Proof. Let £{n) be the index of the largest nonzero integer in n. For every m ^ 1 
we have 

E l[a{k,n,)^ E n«(^'"'«) = n(E«(^'"))- (4.7) 

iiGAA:^(n) ^ m fe=l ni ,...,n„ ^ fe=l fe=l n^O 

The left hand side and the right hand side are clearly increasing in m, and we 
obtain the lemma by letting rn — > cx3. D 



Lemma 
be finite, 
here. 



4.1 



also holds for complex a{k, n) under the assumption J2k > o Sn > i k(^i '^)l 



i'liis can be proved using dominated convergence, but it is not needed 
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We introduce a Ricmann approximation to the critical density (2.9 1 which will 
be useful in Proposition |4.3| below. 



P, 



(A) 



-E«-mE 



3> 1 



|A| 



^-i^{k) 



fc^O 



Lemma 4.2. lim 



|A|- 



(A) 
> Pc == Pc 



(4.8) 



-je(fc) 



Proof. Since e ^'■'^' is uniformly continuous, the Riemann sum |A| ^ X^fc^to ^ 
converges to / e^-''^*^'^) dfc for each j. We need to show that the limit |A| — > oo can 
be interchanged with the sum over j and we use dominated convergence. We sum 
separately over ||fc|j ^ 1 and \\k\\ > 1. For ||A:|| > 1 we use e{k) > c > 0, so that 



M 2^ ^ ^ ^ lAI Z^ 



~e{k) 



|A| 



||fc||>l 



|A| 



(4.9) 



feSA* 



for j ^ 2. The latter sum is easily seen to converge using Eq. (3.5 1, so the right hand 
side is bounded by Ce~^'^/^ , which is summable. For ||fc|j ^ 1, we use e{k) > a|jfc|j'' 
with a > 0. Since e"-'""'^"'' is decreasing, we can estimate it with integrals, namely 



Al 2^ ^ lAI Z^ 



~aj\\kr 



|A| 



k^O,\\k\\ ^ 1 



|A| 



fc/O 



^2° 



2''d7r'^/2r(^) 



(4.10) 



r{l + l)r^{ajY/^' 



The only relevant term in the upper bound is j '^''', which makes the sum over 

ice 

D 



j summable in Eq. (4.8 1. The claim follows from the dominated convergence 
theorem. 



Let A/a = Un ^ o-^a.n be the set of occupation numbers on A* where ^j, Uk is 
an arbitrary finite number. Using Lemma |4.1| we find 



Z=Y^ Y{N)= ^ n e-^^'^"'=V 

JV ^ neAfA fcsA* 

"0=0 



'je(k) 



fe^O n ^ 



J" > 1 



k^O 



(4.11) 



The last identity follows from (4.31 since mik^oe{k) — c(A) > 0. Z is finite by 
Lemma 4.2 This allows to introduce the following probability measure on [0,oo): 



MA 



I E ^(^)^ 



N/\A\- 



(4.12) 



N > 



The motivation for /^a is that the distribution of the occupation of the zero mode 
can be expressed as 



N 



Prob(^^a)= Y. ^. 



j=[aN] 
n(l-a)p 



YjN ~ j) 

YiN) 



L(l-a)WJ 

E 

3=0 



IjV-j 



Yjj) 
Y{N) 



(4.13) 



/;^-"^^fa(|A|(p-.))dMA(.) 

J^'h{\A\{p-s))d^,A{s) ■ 



SPATIAL RANDOM PERMUTATIONS AND POISSON-DIRICHLET LAW 



11 



Here h{x) can be any function interpolating the values ft.„ at n G N, e.g. linear 
interpolation. We use the notation purely for convenience and will never evaluate 
h{x) at non-integer points. 

We now have all the elements that allow to state and to prove the key properties 
leading to the macroscopic occupation of the zero Fourier mode. 

Proposition 4.3. 

(a) /lA ^^ Sp^ weakly as |A| — >■ cx). 



(b) Let A(A) e R such that |A(A)| < |A| 



l-v/d 



,A(A)(X-p<'^)) ■ 



then 
^ 1. 



The parameter rj in the claim (b) is the one that appears in the condition for 
e(fc), see the paragraph after Eq. (2.6). The relevant aspect of the claim (b) is that 
the expectation is bounded uniformly in the domain even though A(A) diverges. 
Markov's inequality then gives the following concentration property for |A| large 
enough, which will be used later: 



l-v/d 



Prob^J|X-pe| >e)< e-'l^l ' S^^(e 



l-v/d 



, 1-v/d 



|A|^^(X-p(*)) _^ ^-IAI^^(X-pW) 



+ e 



<3e-^l^l" 



(4.14) 



Proof of Proposition \4-3\ (a) follows from (b), see (4.14). For (b) we note that 
a{k,n) = e^(^('^)^'^('^)/l^l)" /i„ fulfils the assumptions of Lemma 



4.1 



and thus 



i?MA(e'(^)^) = / eM^)M/.A(s) = 4 E ^^^)^ 

Jo ^ ^^ ^ n 



A(A)iV/|A| 



N > 



y^ TT e"(^('=)~^('^)/l'^l)"'' hn. 



(4.15) 



tisWa fe/0 
"0=0 



Since inik^oe{k) — A(A)/|A| > 0, ( |4.3[ ) applies, and together with (4.111 we obtain 
i?^Je^(^)^) =exp(^ ie-"^ ^ e^^'^W (eJ^(A)/|A| _ i)) (4.I6) 



i> 1 



k^O 



By ( |4.8| and rearranging, we get 

j> 1 



^^^^(g.t.n— Po--')) ^exp(^ Ig-a, ^e~,e(fc)^g,A(A)/|A| _i_j-A(A)^y 



^MA I 



(4.17) 
We show that the exponent vanishes as |A| — > 00 using dominated convergence. 
We use I e^ — I — x\ ^ ^x^ e'^' (which is easy to check using Taylor series) with 
X = jA(A)/|A|. The exponent in the right hand side of (4.17) is bounded, in 
absolute value, by 



1 






3A'(A) ^.7(^e(fc)-|A(A)|/|A|) 
|A| e 



(4.18) 



Since e(fc) ^ a\A\-'^/'^ and |A(A)|2 < |A|i-''/'', the bracket is bounded above uni- 
formly in j for all |A| large enough. The sum over k has been estimated in (4.9) 
and (4.10). Since ^- e~°'' j^'^^^ < cx), we can interchange the limit |A| — )• 00 
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and the sum over j by dominated convergence. The bracket in (4.18) tends to 
as |A| -^ oo, for all fixed j. It follows that (4.18) converges to zero and (4.17) 
converges to one. D 

Proposition 4.4. Suppose that p > p^. Then, in the thermodynamic limit |A|, N — !■ 
oo, 

'l z/a<l-^, 
z/a>l-^. 

Proof. We show that, for every e > 0, 



Prob(^ > a) 



Prob(|^-£^| >e) ^0. (4.19) 



Using the expression (4.12 1 which involves the measure /iA, we write the probability 
as 

J- + J+ 
J-+Jo + J+'' 



Prob(|t-g^|>6)^ ^ r; \ ^ (4-20) 



with 



p-^ h{\A\{p~s)) 

J- = / TJYTT, ^d/iA(s , 

Jo H\Mip~Pc)) 

p+-^ h{\A\{p^s)) 
Jp^~ep ^(|A|(p-pc)) 

J+ = / -TTYvr, ^d^iA s . 



By (4.2 1, the ratios of h{-) are bounded above and below for all < s < p — c 
uniformly in A. This shows that J_ — >■ and that Jq is bounded away from as 
|A| — ^ 00. As for J+, we need to be careful with the integration over s close to p. 
By (4.1 ), the ratio in J_|_ is bounded by 

Zu'"l ^ ^^(^|A|P)^^ (4.22) 

which is valid provided |A|e ^ 1. Then 

J+ < e'^{C\A\p)^^E^,{Uyp^+,p), (4.23) 



and it goes to zero by Eq. (4.14). D 

4.2. No macroscopic occupation below the critical density. For p < pc, the 

support of the Dirac measure to which the measures pA converge lies outside of the 
interval [0, p], and the above argument fails. In its place, we use the formula 

Prob(no ^ j) = -1- V TT e— (*=) h,,, = ^^Zj''\ (4.24) 



' V 




Y{N) ^ 
"0 ^ i 


nSA/A,, 





where 

" '' ' " ■■' (4.25) 
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We apply summation by parts E{f{X)) = /(O) + E^ti [fU) - fU - 1)] P{X ^ j) 
to the function e"'*'""/^ , and we find 



E{e- 



-Xna/N - 



(1 



^\/N 



N 



= 1 



Y{N) 



J2e-^^f^Y{N^j,j) 



3 = 1 
-A /^ X//V\ W-1 



(4.26) 



3=0 



Y{N) 

For the last equality, we used a change of summation index j —> N — j . 
Proposition 4.5. Suppose that p ^ Pc- In the thermodynamic limit \K\, N — >■ cx), 

Prob(^ >(5) ^ 0. 
for all (5 > 0. 



Proof. By (|4.26|), 

''(l-e)Af 



N-l 



N 



3=0 



Y{N) 



E 



y(j,^-i ) 

r(7V) 



(4.27) 



From (4.24| it is obvious that y(N] ^ ^^ ^^ ^* ^^ equal to Prob(rio ^ N ~ j). 
Thus the second term above, along with the prefactor 1/iV, is bounded by conste. 
For the first term, we use the inequality 



h. 



sup ■ 



r+j 



^Coil+jf^ 



(4.28) 



for some constant Co; this follows from (|4Aj) and (|4^. Then, Y{N,j) s$ y(7V)Co(l- 
j)^'', and 






j=0 



Now 



Y{N) 
Y{j) 



3=0 



Y{N) 



-|A|(9A(j/|A|)-gA(p)) 



(4.29) 



(4.30) 



Y{N) 

where qh{p) = ^rST lc)gy(|A|p) is the finite volume free energy associated with 
the partition function Y . It was shown in [4 , under conditions on the coefficients 
aj that are more general than the present ones, that ^a converges uniformly on 
compact intervals to a convex function q, and that p i— > q{p) is strictly decreasing 
for p < pc- Thus for each e > there is 6^ > such that q[^{j/\A\) — q{^{p) > b^ for 
all |A| large enough, and all J ^ (1 — e)N. So 

L(i-«)A'J ^f . .r .. L(i-«)JVJ 

E yl/^ ^goe-"-!^! ^ iN-j + ir^Coe-''^^^''N^-+\ 

3=0 ^ ' 3=0 

(4.31) 
which converges to zero as A^ — >■ oo. Since e was arbitrary, we have shown that 
^( Q>^no/N ■) _^ I fQj. .^YY \>Q^ which imphes the claim. D 
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4.3. Occupation of nonzero modes. We now turn to the modes k y^ 0. Recall 
the bound C{s) for the ration of /i„ in Eq. (4.2 1. 

Lemma 4.6. For < a < I, k E A*, and j ^ 0, we have 

Prob(nfc > j) ^ C{^f e-^'(i-'")^('=) . 



Proof. Using (4.2 1, we have 

ProbK^j) = ^ E n e— ^('^'^ /.„,, 

^ ' tiGA/a.jv fc'GA* 

~YiN) ,^ /i„, li "'=' (4.32) 

g(l/a) ,(i_^)^(fc) Y^ -pr „^,e(fe') 

^ y(iV) Z^ li *^ '^"'='- 

We indeed estimated ''"''^'^~'" ^ C'(^) ^ C'Cl/o")- Since e(/c) > e(0) = 0, we 
get an upper bound by replacing the constraint Uk ^ crj by uq ^ aj. Then 

no ^ ffj 



^W „.^ .. /^"O 



E ^^'^'^^^^ n e-"-^('=')/z„„ 



neA^A.jv ° fc'eA* 

"0 > j 



(4.33) 



Y{N) 
<C(l/(7)2e-^'(i-^)^('^'). 



^ ' tiSA/a.jv fc'SA' 

"0 > j 



D 



We now define three sets of occupation numbers, each of which will be shown 
to have measure close to one. Let i> = max(0, 1 — — ); we will prove in the next 
section that v = v, but we do not know this yet. The sets are 

A,^ [ucNk.n : |^-/>| <e} 

B,,s = {n e Nk,n : Y nk< eN] 

o<j|fci|<5 (434) 

Cf^^s.M = |n e 7VA,Ar : ^ TLk < eN\. 

keA-* ,\\k\\ ^ (5 
rik > M 

Proposition 4.7. For any p > 0, we have in the thermodynamic limit N, |A| -^ cx3.- 

(a) For any e > 0, Prob(Ae) -> 1. 

(b) For any e > 0, there exists S > such that liminf Prob(_Be .5) > 1 — e. 

(c) For any e,S > there exists M > such that liminf Prob(Ce ,5.m) > 1 — e. 
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Proof. The claim (a) immediately follows from Propositions 4.4 and 4.5 For (b), 
we use Lemma [46| with cr = 1/2 to get 



E{n,) = ^ ProbK ^^)^ C[2f ^ e^^C^'^/^ = (^)f/_ (4.35) 

i ^ 1 i ^ 1 

For every (5 > we get, by Markov's inequality, 
Prob(i?.%) ^ ^ E ..vl ^- ^^ / ^^n (4.36) 

0<l|fe||<5 '^ •'ll'-ll<'' 

By the assumption e{k) ^ a||fc||'' with rj < d, the integral is finite, and thus 5 can 
be chosen so small that limsupProb(i3J^) < e. 

For (c), we define F{n) = EfceAMlfe|l ^ 5 ^^fclfnt ^ m}- Note that Prob(C^ ,, j^J = 
Prob(i^ ^ eA^). Now 

E{F/N) = ^ J2 E{nkln,^M) 

feeAMlfcll ^ s , ^ 

1 ' , (4.37) 

^N ^ (MProb(nfe > Af ) + ^ Prob(nfc ^ j)j , 
fceAM|fc|| ^ s j>M 

where the last equality is summation by parts. By Lemma |4.6[ 

EProbK^j)<C(2)2 Y^ e-MW2^C(2)2e-^'^-W/^ ^^(J^_^ . 

j>A/ j=M+l 

(4.38) 
Define c{6) = inf|H;|[ ^ s^{k)- Note that c((5) > for all 6 > follows from the 
properties of e(fc) stated below (2.6). Then, 

^-Me{k)/i (^ j_ 

M(e'^«/2 -1) 

(4.39) 
The sum above, along with the factor 1/A^, converges to a Riemann integral which is 
finite thanks to our conditions on e(fc). Therefore lim sup i?(F/iV) ^ CiM e~'^'^^^^/'^ , 
and Markov's inequality implies that limsupProb(C^|5 j^j) ^ K q~Mc(S)/4, ^ Choos- 
ing M large enough for given e, 5 proves the claim. D 

5. Cycle lengths of spatial permutations 



E{F/N)i^C{2fMe~^'^^'^/^- J2 e-^''^'^^/" (l + ^ 



We now prove the claims of Theorems |2.1| and |2.2[ starting with the fraction v of 
points in infinite cycles. We denote by Prob„(7r) the probability of a permutation 
TT G iSn in the nonspatial model with cycle weights. That is, 

Prob„(^) = T^^l[ e-"^'^^('^) (5.1) 

J ^ 1 



with h„ the normalization defined in (3.9). We also write En for the corresponding 
expectation. We keep the notation Prob, E for probability and expectation with 
respect to the spatial model. 

Recall that we defined D = max(0, 1 — — ). 

Proposition 5.1. Under the assumptions of Theorems\2.1\ or\2.^ we have v — v . 
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Proof. We use the Fourier modes decomposition of Section [SJ Recall that tt = (tt^), 
rjk — rjliTk), and Vj = ^^ r^fc. We have, for fixed K > 1, 

i-H^)>K j>K 

= ^QT.Jr,o)+E{j^ J: E^-^.^+K^ E E^-^.-)- (5-2) 

j>^ 0<||fc||<<5j>-R' \\k\\^Sj>K 

The first term of the right-hand side is equal to 

^(^ E J'^^") = E ^Prob(no = n)i?„(^ E ^^^)- (5-3) 

j>Ar n ^ j>K 



It follows from Proposition 4.7 (a) that ^ — ^ ;> in probability as |A|, A^ — ?► oo. In 
addition, we have 

En[^ E •?''j) = Prob„(4 > if), (5.4) 

where £i is the length of the cycle that contains the index 1. It was shown in [T6jl6] 
that the latter converges to 1 as n — >■ cx). We have thus proved that, for any finite 

lim ^fi- y jr,o) =£^. (5.5) 

i>K 



The second term in the right-hand side of (5.2) is less than E{j^ X]o<||fc||<(5 "fc) ^^^ 
this is as small as we want by choosing 5 small, see Proposition 4.7 (b). The last 
term is less than 

For any (5 > 0, this can be made small by choosing K large, see Proposition 4.7 (c). 
This shows that both vk and vk converge to v as K ^ oo. D 

The next step is to prove that the distribution of cycle lengths of nonspatial 
weighted random permutations is asymptotically equal to Poisson-Dirichlet. 



Proposition 5.2. Assume that aj -^ a as in Theorem 2.1 Then, under the 
probability measure (5.1), 

,»(i) /'(™)\ 

( — ,...,^)^PD(c-"). 

Proof. Let us order the cycles of a permutation tt according to some rule, such as 
their smallest element. That is, the first cycle is the one that contains the index 1; 
the second cycle is the one that contains the smallest element that is not already in 
the first cycle; and so on... Let ^1,^2, ■■ • be the cycle lengths with respect to this 
order. We prove that 



Vn'n — ^1' ^ n — ii — ■ ■ ■ — £m-i - 

converges (in distribution) to i.i.d. beta random variables with parameters (1, e~" ). 
It then immediately follows that (^, . . . , ^) converges to GEM(e~" ), and that 

{- — , . . . , - — ) converges to PD(e^" ). We proceed by induction on m. The case 
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TO = 1 is just the law for — , whose convergence to the beta random variable was 
proved in |161 15] . For ?77, ^ 1 , let 

A - {(ci, . . . , c™) e {1, . . . , nr : t ^ ai, . . . , ^..^J:?.,^,^ ^^ a™}- (5.6) 
Then 

Prob„((^,..., „_,^"J+L,„ ) e^x [0,a„,+i]; 

= ^ Prob„ \J.i= ci,...,(.m^ 
(ci,...,c^)eA 

Prob«( „_£|^"+L<;^ SS flm+l 1^1 = Ci,. ..,4n = C„ J . (5.7) 

It is not hard to check that 

Prob„ f „_^^"+L£^ s^ Om+i I ^1 = ci, . . . , £m = c„ 

= Prob„_,,_..._e„ ( „_,^!^._,^ < a™+i) . (5.8) 
In addition, any element (ci, . . . , €,„) in ^ satisfies 



n — ci 



c„, ^7i]^(l-a,;). (5.9) 



It follows that (5.8 1 converges in probability to the beta measure of [0, a„j+i] uni- 
formly in ci, . . . ,Cm. Then (5.7| converges to a product of beta measures of the set 
x™r[^[0,ai] by the induction hypothesis. D 

Finally, we relate the distribution of long cycles of the spatial model to that of 
nonspatial random permutations. Let 

A = [ai, 61] X • • • X [a™, b„,] C (0, 1)". (5.10) 

Proposition 5.3. If v > Q, we have for any m ^ 1, 

lim Probff ,..., ) eA]= lim Prob„(( ,..., ) eA 

Proof. Let us write i], for the length of the longest cycle of the permutation tt^ 
corresponding to fc € A* . We clearly have 

Probf sup ^ > e) < Probf sup '^ > eV (5.11) 



It follows from Proposition |4.7| (b) and (c) that the right-hand side vanishes in 
the limit |A|,A^ -^ 00. The zero Fourier mode is consequently the only one that 
matters, i.e. 

hm Probff^,...,^) eA) = lim Probf ( \-, . . . , \- ) G A 

= lim Probff^,..., ^^"j gaV (5.12) 



The last identity follows from Proposition |4.7| (a). Since no — > 00 as |A|, A^ — >■ cx), 
the last term converges to the asymptotic joint probability of the m largest cycles 
in nonspatial random permutations with cycle weights. D 
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Theorem |2.1| clearly fo llow s fro m P ropositions |5.1[ |5.2[ and |5.3[ Theorem |2.2 
follows from Propositions 



5.1 



and 



5.3 



and from the fact that — =^ 1 for random 



permutations with cycle weights of the form e "^ = j ^ with 7 > 0, see [B]. Notice 



that Proposition 5.3 is trivial here for m ^ 2, as both sides of the identity converge 
to zero. 
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Addendum 

As Antal Jarai has pointed out to us, the proofs given in this article are not 
complete: We use the Poisson summation formula without assuming decay proper- 
ties of the function; and our assumptions are not sufficient to ensure convergence 
of the Riemann approximations to integrals. 
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We actually need to make two additional assumptions on the function c^^^^' . 
These assumptions seem to hold for all reasonable choices of ^, and in fact we 
have not been able to find an example of a function e~^ that fuUfills our original 
assumptions (in particular to be positive with positive Fourier transform), but not 
the new ones. The main statements of the article stand unchanged. 

We first discuss the Poisson formula. Rather than assuming extra decay prop- 
erty, we use positivity of the function and its Fourier transform. 



First additional assumption. Suppose that 

z&L'i 

is a continuous function of x, for all L large enough. It follows that (c^^'^ )*^ is also 
continuous and integrable (Young's inequality). Next, let us recall the following 
Fourier inversion formula, that is valid for any continuous function / on [—1,1]'^: 

!{x) = lim ^ /(fc) c-'^'^l'^l' e^'^"^-^ , 

for all X in the torus, and where 

Kk) = / /(:c) c-'"'"-" dx. 



This gives, after rescaling, 

(e-«'^)*-'(a;) = ^lim V 



lAU" 



Q-ieC^) Q-^';|fc| g27rifc-x 



For X = we can use the monotone convergence theorem to take the limit r/ — > 0, 
and wc find 



(e-«M*^(o)^E(°"*)*'(^^) = m^ 



c-J£(*=) 



This allows to get Eq. (3.6), from Eq. (3.3). 

Next, in the proof of Lemma 4.2, we mistakenly assumed that Ricmann sums 
of uniformly continuous functions always converge to integrals. There are indeed 
counter-examples, though such functions may not be positive with a positive Fourier 
transform. 

Second additional assumption. Suppose that 



1™ I A I 

L^oo |A| 






It follows from continuity of c '^^'^^ that 

lim ^ V e-^^C"^) = / e-^^^'^M/c 

for all j ^ 1, and monotone convergence allows to take i? — ?■ oo. 
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Finally, our proof that (4.18) tends to uses that e~2'^('=) is Riemann integrable, 
which we have not assumed. But it is not hard to repair this. First, the term j = 1 
in the exponent of (4.17) is easily seen to go to 0. The sum over j ^ 2 can be 
bounded by (4.18) with j ^ 2. The bracket is less than a constant. We get the 
bound 

which is finite, see the proof of Lemma 4.2. Notice that the constant C above also 
depend on the weights e^"^ , but these are bounded. 

We are very grateful to Antal Jarai for pointing out the gaps in the proofs. 



